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Mean field theory for the time evolution of quantum meson fields is studied in terms 
of the functional Schrodinger picture with a time-dependent Gaussian variational wave 
functional. We first show that the equations of motion for the variational wavefunc- 
tional can be rewritten in a compact form similar to the Hartree-Bogoliubov equations 
in quantum many-body theory and this result is used to recover the covariance of the 
theory. We then apply this method to the O(N) model and present analytic solutions 
of the mean field evolution equations for an N-component scalar field. These solutions 
correspond to quantum rotations in isospin space and represent generalizations of the 
classical solutions obtained earlier by Anselm and Ryskin. As compared to classical 
solutions new effects arise because of the coupling between the average value of the 
field and its quantum fluctuations. We show how to generalize these solutions to the 
case of mean field dynamics at finite temperature. The relevance of these solutions 
for the observation of a coherent collective state or a disoriented chiral condensate in 
ultra-relativistic nuclear collisions is discussed. 
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1. Introduction 



Several studies have been devoted to the determination of the evolution of a self-interacting 
scalar field both at the classical and quantum level. Some calculations have considered 
the possible formation of a disoriented chiral condensate in ultra-relativistic nuclear colli- 
sions 1 )" 6 ). The evolution of a scalar field is also important to describe some aspects of 
inflationary models 7 ) " . 

Most analytical results in this domain concern the classical equations of motion. This is 
the case in particular for the rotating solutions in isospin space which were investigated some 
time ago for the non linear sigma model by Enikova, Karloukovski and Velchev 13 ) , Blaizot 
and Krzywicki ^ and Anselm and Ryskin 5 ) . These last authors discussed the observability of 
a classical coherent pion field created during an ultra relativistic nuclear collision in a finite 
volume. They examined the effect it would have on the outside volume where it would act 
as a source. Anselm and Ryskin concluded that a large number of pions, of the order of fifty 
in typical collisions, would be emitted in an interaction vertex, all of them with almost the 
same momentum. Furthermore some events would contain mainly neutral pions while some 
other events would contain in contrast mainly charged pions. 

Classical equations are of interest when the pion density is high enough to make a de- 
scription by a coherent state plausible. An interesting discussion of this limit and several 
related questions can be found in the section on waves and particles in the book by Peierls 
on "Surprises in Theoretical Physics" 1A \ In this section Peierls points out that physicists 
originally encountered only the wave aspects of light and only the corpuscular aspects of 
particles with mass. He argues that this is no accident at all because the analogy between 
light and matter has very severe limitations. In the case of neutral pions for instance, the 
classical field description implies superpositions of states with different numbers i.e. corre- 
sponding to substantial energy changes. This implies rapid oscillations in time in the case of 
free fields. Peierls therefore concludes that the question of the knowledge of the field phase, 
which is central for a classical description, becomes of academic interest in this particular 
case. 

In order to have some deeper understanding on the domain of validity of the classical 
field picture and the coherent or collective effects it implies, it is highly desirable to have 
calculational tools incorporating quantum effects as well as a consistent treatment of in- 
teractions. The most natural framework for this purpose is the well known mean field or 
Hartree-Bogoliubov approach, which assumes that at each time the density matrix of the 
system is given by the exponential of a quadratic form in the field operators (or that the state 
of the system is described by a Gaussian functional in the special case of zero temperature). 
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Several numerical investigations of the mean field equations, which generally involve 
large amounts of computational work, are now available for a wide variety of geometries and 
initial conditions. Some analytical results have also been obtained concerning the return 
of the system to equilibrium in linear response theory 10 ) and the behaviour of mean field 
solutions in a boost invariant geometry at large proper times 15 - 1 . In our recent work 16 ) we 
have used a covariant form of the mean field equations to construct a quantum generalization 
of the time-dependent running-wave-type solutions of Anselm and Ryskin for the classical 
equation of motion of non-linear sigma model. 

The purpose of the present paper is three-fold : i) First we make a systematic review 
of the properties of the equations of motion in the Gaussian variational approach. By 
introducing a generalized (one-body) density matrix, we show that the variational equations 
of motion can be written in the form of Hartree-Bogoliubov like equations, as encountered 
in many-body theory. An attractive property of the density matrix is that its eigenvectors 
are found to be identical to the so-called mode functions. These define at each time a set 
of creation and annihilation operators having the Gaussian state as a vacuum. The mode 
functions satisfy a set of Klein-Gordon type equation in the presence of a mean field and 
they can be exploited to recover the covariance of the theory, ii) The second purpose of this 
paper is to apply the method to the O(N) sigma model and describe a set of time-dependent 
solutions which correspond to rotations in the internal symmetry space of the quantum 
fields. This generalization of the classical Anselm-Ryskin solutions was presented in our 
recent work 16 ) also using the covariant formalism, iii) In this paper we further show (and 
this is the main result of this present study) that the previous solutions can be extended to 
finite temperature. These solutions allow one to investigate the phase diagram of the model 
as a function of thermal excitation energy and/or collective excitation energy. 

The present article is organized as follows. We describe the framework of the Gaussian 
approximation in section 2 in the case of a single component scalar meson field. The equiva- 
lence of this approximation with the time dependent Hartree-Bogoliubov equations is briefly 
reviewed in section 3. In section 4 a covariant form of the mean field equations is presented. 
It uses the so called mode functions, which are the eigenvectors of the generalized density 
operator. These functions have remarkable properties needed later on to extend the Anselm 
Ryskin solutions. In section 5 we consider the case of a scalar field with N components and 
an O(N) symmetry. We review some of the results obtained earlier in this model for the 
classical equations of motion and describe the structure of the mean field equations. We 
show that rotating solutions in isospin space of the mean field equations can be constructed. 
In section 6 an extension to describe mean field equations at finite temperature in the O(N) 
model is given and particular analytic solutions of these equations are also presented. They 
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correspond to quantum rotations at finite temperature in isospin space. They include a self 
consistent treatment of quantum fluctuations around the classical solutions of Anselm and 
Ryskin in the special case of zero momentum and zero temperature. The detailed derivation 
is included in Appendices A ~ D. Numerical results are presented in section 7. We discuss 
in section 8 the importance of quantum effects, with special focus on the self consistency 
conditions occurring in our equations. The relevance of our results regarding the possible 
formation of a coherent collective state or a disoriented chiral condensate in ultra relativistic 
collisions is sketched in section 9. 



§2. Mean field approximation in the functional Schrodinger picture 

We consider in this section a single component scalar field ip described by the Hamiltonian 
density: 

= ^ 2 (x) + \ (Vy?(x)) 2 + (x) + A^( x ) (24) 

where the operator vr(x) is the functional derivative —i8/8ip(x.). The mean field equations 
are obtained by assuming that at each time t the wave functional of the system is a Gaussian 

$ fc(x)] = M exp (V|v? - <p) - {<p - (p\{^ + iE)\<p - , (2-2) 

where G, E, <p, n define respectively the real and imaginary part of the kernel of the Gaussian 
and its average position and momentum. We have used the short hand notation 

(n\ip) = J 7f(x,t)v?(x)dx. (2-3) 

The resulting equations are found to be 

G = 2{GE + EG), 

E = lG- 2 -2E 2 -\(-A + ml + ^ 2 + \G{^)), 

(24) 

= -7T, 

TT = (-Z\ + m 2 1 + ^ 2 + fG'(x,x))v?. 
The vacuum state corresponds to the static solution 

27 = 0, tt = 0, y? = y? , G = G , (2-5) 
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where 

The quantity /j, has to satisfy the so-called gap equation 25 ) 

" 2 = m » + ^ <x| 27^+7? |x> + ^»' (27) 

In the symmetric phase the expectation value of the field ipo vanishes while in the symmetry 
broken phase it must be such that 

m « + I <x| 27=3+7? |x> + ^° = ' (2 ' 8) 

This last equation implies that fi 2 =\(p 2 ] /3. 

The previous equation requires a regularization scheme such as a discretization of the 
Laplacian operator on a lattice with a mesh size Ax = 1/A or a cutoff A in momentum 
space. To make the gap equation finite when the scale A goes to oo a popular prescription 
is to send the bare coupling constant to zero according to the formula 

= T + T^M-)- (2-9) 



2X R A 16vr 2 ° e/i' 

An interesting consequence of this prescription is that it also makes the dynamical mean 
field equations finite 17 ^' n )' 12 \ Note however that one difficulty arises because the bare 
coupling constant becomes infinitesimally small by negative values when the momentum 
cutoff goes to infinity, so that the model is unstable i.e. not really a viable field theory. This 
prescription and its limitations have been extensively discussed in the literature 8 )' 25 \ In 
what follows we will mainly consider, unless otherwise specified, our interacting scalar field 
model as an effective low energy cutoff theory 3 )' 27 \ The model is however a useful one 
which contains attractive aspects. It has in particular a single bound state whose influence 
on meson distributions gives rise to interesting effects 17 \ This bound state occurs because 
the \(p A coupling is a contact force which can accommodate one bound state only 28 ). 

We have in this section considered only pure states i.e. a zero temperature theory. The 
generalization to statistical mixtures and non zero temperatures is straightforward and the 
corresponding mean field equations can be found in the literature 8 ) . 

§3. Dynamics in Hartree-Bogoliubov form 

It is sometimes convenient to rewrite the mean field evolution equations in a form which 
exhibits in a more transparent way the underlying Lorentz invariance of the equations and 
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which furthermore reduces to a set of linear equations for a non interacting meson gas. Such a 
form is provided by the familiar time dependent Hartree- Bogoliubov equations 26 - ) ' 8 - ) ' 29 - ) which 
are particularly suited to investigate boost invariant configurations. In order to construct 
these equations we introduce the generalized density matrix Ai which adequately implements 
the requirements of Lorentz invariance. It is defined by 

M ^ t) + lJf(f^ .), (3-1) 

2 ^ (Tr(x)Tr(y)) -i{n(x)<p(y)) J 

where = {p — , 7r = 7r — n, lis the unit (Dirac) matrix and expectation values are 
calculated with the Gaussian functional &(t). Explicitely 

, J I -2iGE G \ 

M = \ , . 3-2 

\ ± + AEGE 2iSG J 

From its structure it can be checked that the generalized density matrix satisfies 

M 2 = \. (3-3) 
4 

The eigenvalues of the density matrix are thus ± 1/2. Note that if (u,v) is an eigenvector 
of M. with eigenvalue 1/2 i.e. 



-2iGE G 
± + 4EGE 2iEG 





(3-4) 



then (u*,—v*) is also an eigenvector with eigenvalue -1/2. Eigenvectors are conveniently 
normalized to 

J dx (u*(x)u m (x) + <(x)w m (x)) = ±<J m>n , (3-5) 

with a ± sign for eigenvalues ±1/2. It can also be checked that for eigenvalues +1/2, u and 
v components of eigenvectors are related by 



v x 



(^ + 2^)Mx). (3-6) 



The components u of the eigenvectors are often referred to as mode functions 7 )' 9 ). 

For the particular normalisation condition we have adopted the eigenvectors u n and v n 
provide the following spectral decomposition for the generalized density matrix M. 



Z n>0 



U r . 



where the sum runs over eigenstates with eigenvalues +1/2 only. 



(3-7) 



6 



The interpretation of the mode functions can be made transparent from the following 
observation. For a given Gaussian functional & specified by its kernel G and E the operators 



b n = J dx ^v n (x)(p(x) +M n (x) , (3- 



can be seen to have the state \P as a vacuum state for eigenvalues +1/2. Indeed because of 
the relation (3-6) between the components w's and v's we have 

b n m = 0. (3-9) 

Furthermore because of the normalization condition the operators b and b + can be seen to 
have canonical commutation relations. The mode functions thus define a Hartree Bogoliubov 
transformation preserving commutation relations. 

The evolution of the density matrix is governed by the following equations 

d t G = 2{EG + GE), 

d t {G- l /A + AEGE) = -2(rGE + SGT), 

(3-10) 

d t (EG) = fG- 1 + 2EGE - \rG, 

d t (GE) = IG- 1 + 2EGE - \GR 

In this equation r is the mean field 

r = -Z\ + m 2 + ^ 2 + ^(x,x). (3-11) 

A remarkable property of the generalized density matrix is that its equation of motion can 
be written in Liouville- von Neumann form 

iM = [H,M], (3-12) 

where the generalized Hamiltonian 7i has the particularly simple form 

§4. Mean field equations in covariant form 
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The equation of motion can also be written in terms of the eigenvectors (u, v) of the 
generalized density matrix: 

-d;H° ;)(:)■ 

The normalization condition is preserved by these equations. 

Using the explicit form of the mean field operator r we see that the mode functions u n 
satisfy a set of coupled Klein- Gordon type equations 

[u+ml + hp 1 + ^G(x, x)) u n = 0, (4-2) 

while the component v are given by 

v n = id t u n . (4-3) 

This last equation implies that the normalisation condition can be rewritten in terms of the 
mode function u n only as 

i J dx «(x)<9 t M m (x) - M m (x)<9 t <(x)) = ±<J m>n , (4-4) 

which can be recognized as the familiar normalization of the solutions of the Klein- Gordon 
equation 19 ^. 

For the self consistent vacuum the mode functions are 



u k (x) = l/ v / 2^exp(ik.x)/(27r) 3/2 , (4-5) 

where is the self consistent energy ^k 2 + /^ 2 - 

The evolution equations for mode functions form a closed set. This is because the spectral 
decomposition of the generalized density matrix M. implies 

<x|G(t)|x) = i£Mx,t)| 2 . (4-6) 



n 



The evolution equations can be written in a more compact and manifestly covariant form: 

m 2 (x) =m 2 Q + ^(p 2 (x) + ^(x\S\x), (4-7) 

where x = (x ,Xi,x 2 ,x 3 ) and where S is the Feynman propagator in the presence of an x 
dependent mass 

S = ^ -■ (4-8) 
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Indeed the completeness relation of the eigenvectors of the generalized density 
for each time t, reads 



matrix, valid 



£ 

n>0 



U 



(4-9) 



This relation implies the following ones 



and 



^2u* n (x,t)id t u n (y,t) - J2u* n (x,t)id t u n (y,t) = <J(x-y), (4-10) 

n>0 n<0 



^<(x,tK( y ,t) - 5>nCMK(y,0 = o. (4-n) 

n>0 n<0 



From these relations we find that the quantity 

-i ]T (0(t - *')<(x, t) Mn (y, + 0(*' - t) Mn (x, t)<(y, 0) , (4-12) 

n>0 

is the inverse of the operator d tt + T + ie, as can be checked by acting with this operator on 
the previous quantity. As a result 



(x\S\y) = 0( y x -y )J2n>o u n(^ x o)u n ( y y,yo) 



(4-13) 



and 

(x,z |S|x,z > = (x|G(x )|x). (4-14) 

This version of the mean field equations is particularly suited to deal with calculations of 
the linear response of the meson field to an external source 15 ) . 

At this point let us indicate briefly how the previous equations can be extended to finite 
temperatures. The key is that the eigenvalues of the generalized density matrix are constants 
of the motion. They are equal to ± 1/2 at zero temperature. To have a statistical mixture 
one has to consider eigenvalues which differ from this value. For instance one can consider 
eigenvalues of the form f n + 1/2 where the /'s are thermal occupation numbers of bosons. 

§5. The case of the O(N) model 
The O(N) model corresponds to the following Lagrangian density 

C{x) = \d^ a d^ a - l -m\<f - ^A(0 2 ) 2 , (5-1) 
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where 

2 = 0«0a a = 1,2... N (5-2) 

In the special case of the sigma model N=4 the index a corresponds to the isospin quantum 
number. In what follows we will be refer to it as isospin even when N differs from 4. 
In the classical approximation the evolution equations read 

(n + m 2 + ^ ( f ) 2 (x)j ( f )a ^)=0. (5-3) 

5.1. The classical solutions of Anselm and Ryskin for the sigma model 

Plane wave like solutions of these equations were obtained by Anselm and Ryskin 5 ) who 
focused on the special case of the sigma model N=4 and more precisely its non linear version. 
Their solutions have the following structure 

<f) a (x) = A a cos(Lut-k.r + 6 a ), (5-4) 

with the conditions 

£ Al cos(2£ a ) = ]T Al sin(2# a ) = 0. (5-5) 

a a 

For these solutions the square of the field is a constant ipl in space and time 



The corresponding energy density is 



£= 1 -^(k 2 + u 2 ). (5-7) 



For c<j=k=0 (i.e. in the ground state) the constant (pi is such that 



™o + = 0, (5-8) 
while for non vanishing momenta it must satisfy 

^(pl = -ml-k 2 + LU 2 . (5-9) 

We thus see that time like momenta shift the value of the field outside the chiral radius while 
the opposite is true for space like momenta. 

For a vanishing value of the momentum the previous solutions depend, in the case of 
the sigma model N=4, on six arbitrary parameters (A a ,9 a ,uj with three constraints). For 
A\ + A\ 3> 1A\ one has a state with mainly charged pions while the opposite is true for 
A\ + Al 2Aj. It was argued by Anselm and Ryskin 5 ) that the observation of such events 
may be a signature of the formation of a coherent collective state. 
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5.2. Mean field equations at zero temperature 

In the case of the 0(N) model the kernels G and U of the Gaussian wave functional 
become NxN matrices while the center of the Gaussian (often called condensate) becomes 
a vector with N components O , a— 1,2. . . N. In the mean field approximation the evolution 
of (p is governed by 15 )' 18 ) 



A 



A 



A 



D + m + a V 9 + « tr5, ( X ' X ) 5 ab + -^Sabix, x) 



6 



6 



<^(x) = 0. 



(5-10) 



In the previous equation, S(x,x) is a NxN matrix related to the kernel of the Gaussian 
wave-functional 

S ab (x,x) = (x|G a6 (t)|x) 
and the trace runs over the flavour indices. S(x, y) is the Feynman propagator 



S = 



□ + m 2 (x) — %e 



(5-11) 



where the NxN mass matrix is 



ml + ^<p 2 (x) + ^tTS(x,x)j 5 ab 



(5-12) 



The previous equations (5-10) - (5-12) are non-linear because the motion of the condensate 
involves the mass matrix m 2 (x). Note that the first three terms in equation (5-10) cor- 
respond to the classical approximation considered by Anselm and Ryskin. The next two 
correspond to the contribution of quantum fluctuations whose effect is the object of our 
study. Generalization of these equations to the case finite temperatures have been discussed 
in references 20 >' 21 ' . 



5.3. Quantum generalization of the Anselm Ryskin solutions 

Particular solutions of the previous equations which correspond to rotations in isospin 
space at zero temperature were already presented in reference 16 - ) . For completeness we give 
in the present section a short description of the structure of these solutions. They correspond 
to the following form for the condensate 



tp(x) = U(x)t^^ = exp{i(q ■ x)r y } 



( (f \ 



V : / 
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where = (cu, q) and r y is a generator of rotation in the subspace of flavour 1 and 2: 

/ -i ••• \ 
i ••• 
••• 

v ; ; '•'•■) 

The propagator S is of the form 

S(x,y) = U(x)S^(x,y)U\y), 



with 



and 



ip-(x-y) 



(P^i + «)G» M + <?%) — M 2 + is ' 
It can be seen that the previous expressions solve the mean field equations provided the mass 
matrix M in the propagator satisfies 



Ml 



(5-13) 



To have a closed set, the previous equations must be supplemented by the relation satisfied 
by the condensate (p^ 



-q 2 + m 2 Q + ^ + ^trSM(x,x) + ^(x,x) 



V? = 0. 



(5-14) 



The above solutions depend on the 4 parameters (u, q). For fixed values of these parameters 
the above equations provide the values of the matrix elements M a b and of the condensate 
V?o- The matrix elements M a ^ are found to be diagonal in the isospin index i.e. M a b= M a S a b. 
For rotations in the subspace of flavour 1 and 2 one has M%= M^= \x and one thus needs 
only to determine the four quantities Mi, M2, /i and y?o- In reference 16 -* it was shown, using 
a perturbative analysis, that a disappearance of the condensate <^ , reflecting a restauration 
of chiral symmetry, occurs for u=0 at values of q corresponding to an energy density of 
about (123 MeV) 4 . We refer the reader to reference for a discussion of this result. In the 
next section we focus on an extension of this calculation to the case of rotating solutions in 
isospin space at finite temperature. 
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§6. Mean field equations at finite temperature in the O(N) model 

The mean field equations at finite temperature in the O(N) model are presented and 
discussed in Appendix A. The construction of rotating solutions in isospin space at finite 
temperature is presented in Appendix B. These solutions are constructed in two successive 
steps. The first step is to write the mean field evolution equations in a rotating frame in 
isospin space. The second one is to look for static solutions in the rotating frame. In the 
present section we will just give a short description of the structure of the solutions. As 
compared to the previous section, devoted to the zero temperature case, the basic structure 
of the coupled gap equations (5-13) and the condensate equation (5-14) is unchanged. The 
only difference is that the propagator in the rotating frame now involves boson occupation 
numbers: 

S(p) = (2n(p ) + 1)t V2 7 1 V2 tt^—. , (6-1) 

(PO - WTy) 2 - (p - (\Ty) ~ M 2 +t€ 

where 

= e^T' (6 ' 2) 
as shown in Appendix C, see Eqs.(C-22) and (C-18). Here, uo is the frequency of isospin 
rotation and q is its wave number. The appearance of the occupation number (6-2) breaks 
the Lorentz invariance of our previous results for the vacuum case. The non-zero value 
q means the space-oscillation is realized in isospin rotating frame where the condensate is 
static and uniform. The above expression contains a specific information of the boundary 
conditions by the e-prescription for the solution of a generalized Klein-Gordon equation as 
shown in detail in Appendix D using the mode-expansion of the propagator, (D-3). It is 
shown that the gap equations can be written in the form as in the vacuum (5-13) with 
fluctuation terms involving the matrix: 

j^-^S ab (p;x = 0) , (6-3) 

where 

S(p;x = 0) = - / ^(2n(p ) + 1)7 , 2 , 1 ~ 2 77^-— ■ (6-4) 

J 2m (p — uJTyY — (p — q.T y ) 2 — M 2 + te 

The integral with respect to p can be carried out analytically: it picks up the poles of the 
integrand in the lower complex p -plane as in the vacuum case. We present below the results 
for the pure time-like rotation ( q 2 = uj 2 > , q = 0) and for the pure space-like rotation 
(g 2 = — q 2 < , uj — ). For both cases, off-diagonal matrix elements of S a b(p,x = 0) 
vanish after integration over p or p. This implies that the matrix M 2 b is diagonal as in the 
case for the vacuum. 
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The gap equations can therefore be written explicitly as 

A A N A 
Ml = ml + -v?q + -J2 S cc(x,x) + -S u (x,x) 

c=l 6 

A A N A 
Mj = ml + -<pl + S cc(x, x) + -S 22 (x, x) 

D c=1 6 

= Ml = M 4 2 = • • • = M 2 

A A N A 



while the condensate equation becomes 



(6-5) 



A A ^ A 

g^o = ^ 2 - q 2 - m o - g Xj 5cc ( x ' x ) ~~ g^n^j x ) • ( 6 ' 6 ) 

By using the mass M 2 , the above condensate equation can be recast into a simpler form: 

^jj = -w 2 + q 2 + M?. (6-7) 

The explicit forms of the integrands of the diagonal matrix elements of S a b(p, xq = 0) 
are given below for two specific cases. 

Case 1) Pure time-like isospin rotations (q = 0): 



S u {p,x = 0) = (2n(£+) + l) 



+ 



4uj 2 + M 2 - Ml 



4E + 4E + (E 2 + -E 2 _) 

1 4lu 2 + Ml - Ml' 



+ (2n(E_) + 1) 



4E_ 4E_{E 2 + -E 2 _) 



S 22 (p,x = 0) = (2n(E + ) + l) 



1 4cu' 2 + Ml - M'l 



+ 



4E. 

+ (2n(£_) + 1) 



4E + (E\ — E 2 _\ 

1 4lu 2 + Ml - Ml' 



4E_ 



4E_(E 2 + -E 2 ) 



where 



(6-8) 



■i , M 2 + Ml , 2 



(6-9) 
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We note that E± are also eigenvalues for H(u) in (B-3). Also, the matrix elements S a b(x,x) 
for a, b > 3 are given as 

S ab (x, x) = 5 ab J ^ 3 r^-i ( n ^P 2 + /* 2 ) + I) (for a,b — 3,4, - • • ,N) . (6-10) 



Case 2) Pure space-like isospin rotations (u = 0): 



The expressions of the diagonal matrix elements, Sn(x,x) and S22 ( slightly 
different from those derived in the time-like case, while those of S cc (x,x) for c > 3 are 
unmodified: 



where 



S n (p,x = 0) = (2n(E + ) + l) 



1 M 2 - M 



4E + 4E + (E\-E 2 ) 

1 M 2 - Ml 



S 22 (p,x = 0) = (2n(E + ) + l 



+ (2n(E_) + 1) 
1 



4E_ 4E_(E 2 + -E 2 _) 
M\ - M 2 



4E + AE + (El-E 2 ) 
+ (2n(E_) + 1) 



1 M\ - Ml 



4E_ 4E_[E\-E 2 _) 



^±(p) = p +q + 



2 , .Uf + M| 



± 



(Ml-Mf 



+ 4(p • q) 5 



(6-11) 



(6-12) 



>7. Numerical results 



We can numerically solve the gap equation (6-6) with (6-3), (6-8) and (6-10) for q 2 = 
u 2 > 0, or (6-3), (6-11) and (6-10) for q 2 = — |q| 2 < in a self-consistent manner. We set the 
number of isospin N = 4. A three-dimensional momentum cutoff A — 1 GeV is introduced. 
For this cutoff, the model parameters used here, A and m 2 ,, are taken so as to reproduce the 
chiral condensate ipo to be the pion decay constant f w = 93 MeV and Mi to be the sigma 
meson mass M a = 500 MeV at q^ = and T = 0. Then, condensate equation at T = q^ = 
gives A = 86.7 and m 2 , = — (1019. 5MeV) 2 . One-dimensional numerical integration for |p| 
in Eq.(6-8) is performed for q 2 = oj 2 > case. For q 2 = — |q| 2 < 0, two-dimensional 
integration is required because the square of the inner product (p-q) 2 appears in (6-11) with 
(6-12). However, we replace cos 2 9 with its angle-averaged value {cos 2 9) = 1/3. We safely 
perform one-dimensional integration for |p| with the above-mentioned replacement because 
the resultant curves for condensate (p versus |q| are indistinguishable at T = 0. We checked 
numerically that there is no visible difference at T = 0. 
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Fig. 1. The chiral condensate ipo with q 2 is depicted as a function of temperature. The vertical axis 
represents the chiral condensate tpo (MeV) and the horizontal axis is temperature (MeV). The 
solid curve corresponds to the case = 0. Upside five curves from the solid curve correspond 
to the cases yfcp = lo =20, 40, 60, 80 and 100 MeV, and downside two curves from the solid 
curve correspond to the cases \/—q 2 = |q| = 20 and 40 MeV. 



The chiral order parameter with a certain q 2 is shown as a function of temperature in 
Fig. 1. The increase of temperature leads to the decrease of the amplitude of the chiral order 
parameter <po at each q 2 . Namely, chiral symmetry is partially restored at finite temperature. 
At a certain critical temperature T c , a solution of the gap equation does not exist suddenly. 
At T > T c , only trivial solution tp — exists. The order of phase transition is like the 
first order. In our numerical calculation, T c 's are 97.5, 101.3, 110.9, 123.1, 136.3 and 149.7 
MeV for ^q 1 = uo =0, 20, 40, 60, 80 and 100 MeV, respectively, and 93.4 and 76.9 MeV for 
\J —q 2 = |q| = 20 and 40 MeV, respectively. 

Let us consider the case of time-like q 2 , that is, q 2 = uo 2 > with |q| =0. As uo increases at 
each T < T c , the condensate fo also increases. In this case the condensate rotates uniformly 
with angular frequency uo in isospin space. As a result, the chiral symmetry breaking is 
enhanced. 
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Fig. 2. The chiral condensate tpo with temperature T is depicted as a function of three-momentum 
|q|. The vertical axis represents the chiral condensate (fo (MeV) and the horizontal axis rep- 
resents the magnitude of three momentum |q| (MeV). The solid curve corresponds to the zero 
temperature case. Three curves from the solid curve correspond to the cases T =40, 60 and 80 
MeV. 

In the case of space-like q 2 , that is, q 2 = — |q| 2 < with u = 0, the increasing |q| results 
in the decreasing ip at each T < T c . Further, the critical momentum |q c | exists as is shown 
in Fig. 2. Namely, the gap equation becomes to have no non-trivial solution at a certain 
value |q c |. The transition seems the first order phase transition. For example, |q c |'s are 50.0, 
49.9, 49.7, 47.2 and 37.7 MeV for T = 0, 20, 40, 60 and 80 MeV, respectively. 

It is interesting to compare these critical momenta for space-like condensate to the cor- 
responding values in the classical limit. In classical case, |q c | 2 = Ml/2 is expected at T = 
from the classical gap equation. However, the quantum effect leads to more rapid change 
of chiral condensate as was pointed out by perturbation theory in the previous paper. 16 ) 
This is expected because quantum fluctuation smears out the effective potential and makes 
symmetry breaking more difficult to reach. The coupling to the quantum fluctuation works 
to suppress the appearance of static condensate with longer wavelength. This would have 
important implications for the dynamics of chiral condensate in high energy nuclear colli- 
sions. 
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§8. Importance of quantum effects 

In the work of Anselm and Ryskin, which uses the framework of the classical equations of 
motion, it was pointed out that the existence of rotating solutions in isospin space may have 
observable consequences in ultra-relativistic collisions. For instance it may be possible to 
produce a classical state of the pion field in which all pions nearly carry the same momentum, 
which would manifest itself by the observation of pion "jets" . Since Anselm and Ryskin's 
analysis was restricted to the classical picture, an important question is to find out how 
the previous scenario is affected by quantum fluctuations. A possible way to evaluate the 
importance of these effects is to compare the rotational contributions to the energy density 
arising from classical and quantum terms. The energy density at a given momentum is indeed 
proportional to the pion number density times the pion energy. In the present framework 
the rotational contribution £ R to the energy density (due to non-vanishing E and 7f ) at zero 
temperature is 

S R (x) = iyr 2 + 2 Trace(x\EG£\x). (8-1) 

(See, (A-3) with £ = a = at zero temperature in Appendix A.) The first term is of classical 
origin while the second term is purely quantal. Using the fact that G a b = S a b with n(E) = 
at zero temperature and by solving the equation of the second line in (A- 11) to get E with 
iG^ 1 being iu[G~ l , r y \ in isospin rotating frame, this contribution can be rewritten as 

E R (x) = \uo 2 (y 2 + /) , (8-2) 

with 

j_ f dp (S 22 (p,Xo = 0)-S 11 (p,x = 0)) 2 

J (2vr) 3 S 22 {p,x = 0) + S 11 {p,x = 0) 1 ' 

Let us evaluate this expression in the vicinity of uo ~ at T = 0. In this case we replace the 
kernels S by 1/2^/p 2 + M 2 with the result 

1 * ^ I (I / ? > (84) 



where M\ = M 2 + /i 2 and Mf = M 2 — /i 2 . As a first approximation we replace the square 
masses in the denominator by their averages M 2 so that 

_ = _t A 

(2tt) 3 ( p 2 + M 2 f 

This gives the approximate contribution to the rotational energy density 

1 

,'- i - - 

24ti 2 M 2 J ' 



J~ 1„4 f d P 1 ^ ( o 5) 

- ^ J (2tt) 3 r„2 , M 2^/2 24tt 2 M 2 ' 1 ; 



S R (x) - ^ 2 U + . (8-6) 
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Taking the value of the condensate y?o to be the pion decay constant f w = 93 MeV, M\ 
to be the sigma meson mass M a = 500 MeV and M2 to be of the order of the pion mass 
M w = 140 MeV we find that quantum effects give only a small contribution ( 4.8 %) to 
the energy density. We can also evaluate the quantum contribution (8-3) directly in our 
numerical calculation with three-momentum cutoff. We conclude that the quantum effects 
yield only a very small contribution about 2.8 % in the range of u = ~ 100 MeV. 

We thus expect the number of pions outside of the condensate to be small. This implies 
that, in this case, the discussion of Anselm and Ryskin is unaffected by quantum fluctuations. 
Let us recall however that quantum fluctuations still play an important role since they drive 
the chiral radius away from its vacuum value as excitation energy increases. 

§9. Discussion 

In the present paper we have shown that it is possible to construct analytically quantum 
coherent states which are solutions of the mean field evolution equations in the case of 
an N-component scalar field. We have proceeded in two steps. We have first constructed 
rotating solutions in isospin space in the zero temperature case. In a second step we have 
generalized our results to the case of finite temperature. These solutions represent the 
quantum generalization of the classical solutions obtained earlier for the sigma model by 
Anselm and Ryskin 5 \ 

As compared to the work of Anselm and Ryskin we have found that new effects arise 
as a result of the coupling between the motion of the mean (classical) value of the field 
and its quantum fluctuations. In particular restauration of chiral symmetry is affected by 
quantum fluctuations. Another difference is that, while classical solutions correspond to 
a distribution of mesons with almost the same momentum, a new component containing 
all possible momenta appears in the quantum solution. However we have found that the 
number of mesons outside of the condensate is small. As a result the conclusion of Anselm 
and Ryskin about the observability of the coherent classical state remains unaffected. 

Further questions in the line of the study by Anselm and Ryskin are worth additional 
investigations. For instance it would be interesting to consider an initial configuration in 
which the coherent solution we have built is present in a finite volume. By considering this 
field as a source for the external volume it would be possible to find out the distribution of 
the emitted pions as discussed by Anselm and Ryskin 5 ) . 

Acknowledgements 



19 



One of us (D. V.) wishes to thank Prof. Yosuke Nagaoka and the members of the Yukawa 
Institute for Theoretical Physics of the Kyoto University for the hospitality extended to his 
visit. He also wishes to express his appreciation to the Japan Society for the Promotion 
of Science for support to visit the Yukawa Institute. T. M. is grateful to the members 
of Division de Physique Theorique, Institut de Physique Nucleaire, Orsay, for hospitality 
during his visit. Y. T. would like to express his thanks to the members of the Laboratoire de 
Physique Theorique des Particles Elementaires (LPTPE) of Universite Pierre et Marie Curie 
for the warm hospitality. He also wishes to express his sincere thanks to Nishina Memorial 
Foundation for financial support to visit the LPTPE. 



Appendix A 

Mean field equations of motion at finite temperature 

Let us consider in this Appendix the quantum and thermal generalization of the classical 
equations of the O(N) model. Following the Eboli, Jackiw and Pi prescription, 8 ) let us 
introduce a Gaussian density matrix in the functional Schrodinger picture : 

p[<Pi,<P2] = A/" G exp^-i(7r|v?i - <p 2 ) - (fi - PK^q - i£)\tp\ - <p) 

(A-l) 

where Nq is a normalization factor. Here ( is a mixing parameter and we divide it into 
real and imaginary parts as ( = £ + Aiy/GaVG.The real part, £, is a symmetric kernel 
(£(x, y) = £(y, x)) and the imaginary part a is an antisymmetric one. If ( — 0, the Gaussian 
density matrix is expressed in terms of pure states by p[</?i,</? 2 ] = ^ r [v?2]^'*[v 9 i]- Averaged 
values are easily calculated by 

(O) = Trace(pO) = J J V^V^ p[ Vl , V2 ]0[ V2 , Vl ] , (A-2) 

where 0[(p 2 ,(pi] = <</p 2 1 0| v^i) • 

The energy density calculated for such a Gaussian density matrix is found to be 

(W(x)> = £ (x) + fTrace(x|G'-5(l+£)G'-5|x> 

+iTrace(x| {-A + ml) {G* j^G^) |x) (A-3) 
+ 2 Trace(x|(i: + a)G*±;G*(E - <r)|x) + (V(x)), 



20 



where 

£o = \f* + \{V^ a f + \mlf + lA(v, 2 ) 2 (AT) 

is the classical energy density and the Trace is to be taken over the isospin index a. The last 
term in this equation is expressed as 

(F(x)) = ^A(2y9 2 (^ a (x)y9 a (x)> + (0 o (x)0 o (x))(&(x)&(x)) 

(A-5) 

+4^a(x)(^„(x)^ 6 (x))^ 6 (x) + 2(<£ a (x)£ 6 (x))(<£ a (x)£ a (x))) , 

where (</? a (x)</? 6 (x)) = [G 1 / 2 (l — £) _1 G 1//2 ] a 6(x, x) and repeated indices are supposed to be 
summed over. 

The classical equations of motion for (p a {pc,t) and 7f a (x, t) in quantum and thermal fluc- 
tuations are written as 

^a( X ^) = -*o(x,t), 

(A-6) 

# a (x, t) = (-Z\ + m 2 + f^ 2 + £(0 c (x)&(x))) ^ a + §(£ a (x)£ 6 (x))£> 6 . 

Let us introduce the generalized density matrix in O(N) model as is similar to that of 
4 -theory : 

M ab (x,y;t) = (x\M ab (t)\y) 

(7r a (x)7r 6 (y)) i(7r a (x)(^ 6 (y)) / 2 ' 

where (p = ip — (p , %; = ix — ft, lis the 2N x 2iV unit matrix and expectation values are 
calculated with the Gaussian density matrix p\p\, ^2]- As in the case of zero temperature, 
eigenvectors of this matrix allow one to build mode functions. These functions define new 
creation and annihilation operators via Eq.(3-8) which can be used to bring the density 
matrix (AT) in canonical form. 

In the mean field approximation, we postulate that at each time t the state of the system 
is described by a Gaussian density matrix of the form (AT). Evolution equations in this 
approximation scheme are obtained by considering the Heisenberg evolution equations for 
the operators ipip, ipn and nix and by calculating expectation values in the Gaussian state 
(AT). As a result we find 

iM ah (t) = [H,M(t)] ab , (A-8) 

where we define Ti as 
H = 




21 



/ , X N A N \ 

r ab = \ -A + ml + - ^c(x)^ c (x) + - V(<£ c (x)<£ c (x)) 5 ab 



+ ^ ((^a(x)^(x)) + V?a(x)v? b (x)) . 



(A-9) 



Thus, we can get the equations of motion governing the time-evolution of quantum fluctua- 
tion in thermal fluctuation from the above equation of motion. Noting the explicit expression 
for the generalized density matrix, 

M ah 



-2i [VGj^VG(E - a) 
75(1 + + + ^)VG^VG{E ~ °) 



i i 

4 



ab 



2i [(27 + a)VG-^VG 



ab 

(A-10) 



we can obtain the equations of motion for G, 27, £ and a as 

1 



^2 , ^ 



27 =g(G"^- if) -2(27' 
CT 1 = -2(G- 1 27 + 27CT 1 ) + 2(^ - (777) 
77 = -2(7/27 + 27?7) + 2{G~ 1 a - aG' 1 ) , 
a = -2{aE + 27a) + linG' 1 - G^rj) , 

o 

where we define 77 = G~ l l 2 ^G~ l l 2 . 



(A-ll) 



Appendix B 

Equations of motion in iso spin-rotating frame at finite temperature 



In this Appendix we show how to construct solutions which correspond to rotations of 
the components ipi and ip 2 in isospin space with an angular frequency u. These solutions are 
similar to those developed by Thouless and Valatin 22 ) in the context of mean field theory for 
deformed nuclei and the restoration of rotational symmetry when it is spontaneously broken 
in the mean field ground state at zero temperature. Related work concerning the case of a 
multicomponent scalar field can be found in references 23 ^' 2A \ 

To construct these solutions we first consider the equations of motion in the isospin 
rotating frame. In this frame the Hamiltonian is 

(B-l) 



H Iot = H (N) -Ldab j d xy? a (x)7r 6 (x) , 

where we define oj a b as 0J12 = —0021 = oj and the otherwise are equal to 0. We can also derive 
the equation of motion for Ai in the isospin rotating frame : 

tM ab =[H(uj),M(t)} ab , (B-2) 
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where we define H(u) as 



n(u) = 



-UTy 

r 



-LOT,,, 



-i 

1 




(B-3) 



V J 

and r a b is the same as (A-9). Hereafter, we consider a situation that the time-evolution of G, 
E etc. is only originated from the rotation in the isospin space. In this case, the generalized 
density matrix M. is "static" in the isospin-rotating frame, namely, we obtain 



[H(u),M] = 



(B-4) 



Thus, Ti(ou) and Ai are simultaneously diagonalizable. The equations for Q = G" 1 , E, r\ (£) 
and a are obtained by replacing Q into iuj[Q,T y ] in Eq.(A-ll). For the condensate or mean 
field, we obtain the following equation in the isospin-rotating frame : 



w 2 <Pa 



~ A + m o + \p 2 + ^(^c(x)^c(x)) j 5 ab + ^(0 o (x)&(x)) 1 (p b . (B-5) 



Appendix C 

Spectral decomposition of generalized density matrix ■ 



In this Appendix, we give a spectral decomposition for the generalized density matrix M. 
in the a = case in terms of mode functions, namely, the solution of the Hartree-Bogoliubov 
equations. 

We will start with a situation that the condensate or mean field cp does not depend on 
space and time in the isospin-rotating frame, pointing in the direction of the first isospin : 



ol 



(C-l) 



As is shown later the averaged values (y3 a (x)y3b(x)) have a diagonal form in the isospins 
space. The mass matrix M% b from r in Eq.(A-9) reads 



r, 



ah 



Ml 



A 5 



ab 



M 2 ab , 



/ A N A N \ A 

m l + g S(^c(x)) 2 + - £<^c(x)^c(x))J 5«b + -((^a(x)^ 6 (x)) + 
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= M 2 J ab , 

A A N A 

M a = m l+QVl+Q E(^( X )^W) + 3 «^a(x)^ (x)) + tf&i) . (C-2) 

Further, we will show later that (<£ c (x)<£ c (x)) for c > 3 have same values because isospin 
rotation occurs only in c = 1- and 2-components. Thus, we can set up the masses as 



M* = At* - , 

Ml = Ml = ■ ■ ■ = M 2 N = /i 2 , (C-3) 

where we define Mq and as 

A A A N 

Ml = m 2 + -ifl + -((^i(x)^i(x)) + (^ 2 (x)^ 2 (x))) + - ^(<£ c (x)<£ c (x)) , 

l4 = + ^i(^i(x)) - (v3 2 (x)v3 2 (x)» . (C-4) 
Then, the matrix r has a simple form : 

r ab = {-A + m 2 a + //fo)5 a6 , (C-5) 
where (r z )n = — (r z ) 22 = 1 and the others are equal to 0. Here we define m\ = m\ = M, 



and m\ = fj, 2 for a > 3. Using the above simple expressions, we can easily get the eigenvalues 
for Ti{uS) in momentum representation from the following eigenvalue equations : 

~2 1 ) f U * ) =*(k) f "*) , (C6) 

r(k) -cjt^ J \ v k j y v k j 

where (-Uk, t>k) and E(k) are eigenvectors and eigenvalues for Tt(u) in momentum represen- 
tation, respectively. As a result, we get 



£ 2 (k, a = ±) = k 2 + Ml + u 2 ± yjf4 + 4cu 2 (k 2 + M 2 ) , 

£ 2 (k,a > 3) = k 2 + /i 2 , (C-7) 

where +, (— ) means a — 1, (2). 

In order to get the spectral decomposition for the generalized density matrix M. in terms 
of (wk, Vk) , it should be noted that the matrix M. is expressed with the help of three matrices 

as 

' i./IH o 



.Vt = i;[ 2 VW ^ ) H ^ (C-8) 
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where 



V„ 



G^ 



2i(E + a)G^ + G-^^^- 2t(E + a)G-2 -G--2^l 

i 



\G--2 + 2t^=Gh(E - a) 



Gi 



(C-9) 



We will now look for solutions in which the imaginary part of the mixing parameter, a, 
is equal to 0. This corresponds to the "static" case in the rotating frame. Then, the matrix 
M. is simply expressed as 



M = V 











2V 1"« 



v-\ 



(C-10) 



where V = V o - =0 and V x (= W / (J= o) is an inverse operator matrix of V. Then, for positive 
eigenvalues of M, we obtain 



M 



k(k)>" 
l«c(k)> , 



1 l + &(k) /> c (k)) 
2\|l-£ c (k) \> c (k)) 



il«c(k)> 
il«c(k)> 



= V 



kc(k)> 




, (C-ll) 



where £ c (k) and |iu c (k)) are eigenvalue and eigenvector for the operator £(£), respectively. 
Here, n c are normalization factors. It should be noted that the |u c (k)) is related to |w c (k)) 
as follows : 

k(k)> 



1 



\-?-= + 2iE k(k)) , 



(C-12) 



,2 V / G V \/G 

where |w c (k)) = n c G^ |w c (k)). From (C-ll), we obtain the spectral decomposition for the 
generalized density matrix M. : 



A? 



k,(B>0) c=l 



+ 



k(k))\l l+jc(k) 
| Wc (k)>; 2^1-&(k) 

1 |l + f c (k) 



K(k)) 

-k*(k))y 2^1-Cc(k) 



((^(k)|,K(k)|) 
«(k)|,«(k)|) 



(C-13) 



where we used the relations 

M | 



<(k)) 
>c*(k)) 



1 + &(k) / K(k)> 
2^-^) V-k*(k)) 



(C-14) 



for the eigenvectors of negative eigenvalues. The notation E > means that the sum runs 
over the eigenvectors with positive eigenvalues. It should be noted that H(uj) and M are 
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simultaneously diagonalizable because the relation [H(uj),Ai] = is satisfied in the rotating 
frame. Applying adequate normalization factors n c , we get the ortho-normalization and 
completeness relations as 

J rfx{t» c (k,x)*M c (k',x) +w c (k,x)*t; c (k , ,x)} = ±5 kk , (C-15) 

with + (— ) sign for positive (negative) eigenvalues, and 

,&,?{ ( SS ) ( wk)l ' Wk)i) - ( S» ) ( -<" :(k)i - m }= 1 ■ (c - i6) 

Finally, it is necessary to determine £ c (k) in Ai, while the eigenvectors are directly 
calculable. Since we treat the case a (the imaginary part of the mixing parameter) is equal 
to 0, the generalized density matrix has a simple form in Eq. (C-10). Thus, we can write 
the entropy from this density matrix in the unit ks = 1 as 

S = Y^i 1 + n k, c ) ln(l + n k>c ) - n kiC lnn kiC ] , (C-17) 

k ,c 

where the n kjC 's are defined by 

1 



rc k 



1 + ^ c(k) -1) . (C-18) 



2 \\ 1 - &(k) 

Thus, £ c (k) or n kjC occurring in the generalized density matrix can be obtained by imposing 
that the free energy F(u) in the rotating frame 

F( W ) = <WM>-| (C-19) 

is stationary because all variables are time-independent in the rotating frame. Here, (3 = 1/T 
is inverse of temperature. The variation of (TL(uj)) with respect to n kjC is easily obtained by 
noting that the change 5M. is given by 



5M = J2 ( „ Wk,cMk)> c (k)*). (C-20) 

k, c W k )/ 



Since mode functions (u c (k),v c (k)) are eigenstates of the matrix H(u) due to 
[H(u),M] = 0, we have 

8(H(w)) = Trace(H(u;)5M) = £ £ k , c <5n k , c . (C-21) 

k,c 

Moreover, since we obtain SS = Z) k>c ^ n k, c ln((l + n kjC ) / n kjC ) , thus 5F(uj) = gives us 

n kc = , or £ c (k) = — -4 • (C-22) 

' - 1 ' S V ; cosh/^e v ; 

One can understand that n kjC are nothing but boson distribution functions. 
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Appendix D 

Gap equation in the rotating frame 



In this Appendix, we give solutions for (</? a (x)</? 6 (x)) in terms of the mode functions. In 
a first step we consider the special case of a Gaussian density matrix with zero momentum. 
To treat the most general case of a Gaussian density matrix with finite momentum we 
exploit in a second step the evolution equations for mode functions to generate in a natural 
way solutions with a finite value of the momentum. We can get the results with the finite 
momentum by following the method developed in §5. The gap equation to determine the 
chiral order parameter ip is derived in the closed set of equations. 

From the spectral decomposition of M, we can write (<£ a (x)<£ 6 (x)) explicitly as 

(&(x)ft(x))=A<S*(x,x) 



k,(E>0) * 2 \ 1 -&r(£k,a) 



ak 



x)^ k (x)+ M a * 



ak 



o-k( 



(D-1) 



We can obtain (<£ (x)<£&(x)) by calculating directly the mode functions w" k (x) for TC(u) such 



as 



(6,(x)&(x)> = <u X! Tp^F ( nk,a + D (for a ' b = 3 ' 4 ' ' ' ' ' N) ' (D ' 2) 



where /i 2 is defined by Eqs. (C-2) and (C-3). For a > 3, we can easily obtain the above 
expression by calculating the mode functions directly. Although the isospin-rotating com- 
ponents, (<£ (x)<£&(x)) (a, 6 = 1,2), are easily calculable through mode functions, we give 
another derivation here. The two methods, of course, give identical results. 
Let us introduce the following N x N matrix operator : 



(x,x \S ab \y,y ) = 6(x 

+% - x 



yo) E E \ T^Yck) 



<k( x > z )i4*k(y,yo) 



k,(E>0) a \ 1 ~~ ^( k ) 



U 



crk 



x > ^o)«*k(y> j/o) 



(x|S , a6 (a;o,l/o)|y) , 



(D-3) 



where the mode functions w^ k (x, x ), together with i£ k (x, x ), satisfy the following equations 
of motion in the isospin-rotating frame : 



id. 



« CTk (x, aj J 
w CTk (x,x ) 



— u;t„ 



r(k) -UTy 



u ak (-x,x ) 
v ak (x,x ) 



(D-4) 
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where r has been already defined in Eq.(C-5). Eliminating v^^x, xq), we get the following 
Klein-Gordon-type equations of motion : 

(-d 2 XQ -r + u 2 + 2iuT y d X0 )u ak (yi, x ) = . (D-5) 

The matrix operator S is related to (</? a (x)</? 6 (x)) through the relation 

(x, x \S\x, x ) = M 12 (x, x) (D-6) 

and Eq.(D-l). Here we note that M CT k(x, xq) are also eigenvectors for the operator £(xo) from 
the relation \u c ) oc G 1 / 2 !^) because £(a?o) and G(xo) commute due to the equation of motion 
for the imaginary part of the mixing parameter, a, in the case a — 0. We introduce the 
Feynman propagator 

/ i c I \ _ / i (i+iM v% / \ 

(x,xo|^|y,y ) = (x|^ r -^- y j ^• 2/0 HRwj ly) 
= 0(x o -y ) E <k( x > x )u%(y, y ) 

k,(B>0) CT 

+%o-^o) E E<k(x,^o)^ k (y,yo) , (D-7) 

k,(B>0) <? 

which satisfies 

(d 2 X0 + T-UJ 2 - 2iuTyd Xo )(x,x \S\y,y ) = -i5\x - y) (D-8) 

by using the equations of motion for u^^x., x ) and i> CT k(x, rc ) and the completeness relation. 
Thus, it is symbolically written by 

( x > arol^ly, Vo) = a2 , = ,, 2 ~* a (D-9) 



<9 2 + r - to 2 - 2iujT y d x 



it 



xq M ^^'y^x 

with infinitesimal positive constant e. From the Fourier-transformation for the propagator 
(x.,XQ\S\y,y Q ), namely 

S(x,y) = J ' -^S(P>~ iP{X - y) > (D-10) 



we obtain 



s{p)= MwfoO 



In this way, we get the Fourier transformation of (x, x |5'|y, y ) : 

s(rt = {2nM + 1) ( P „-^-^-^-^ + ie • < D12 > 
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where 



n(p ) = 



1 



(D-13) 



e Ppo - 1 

with the help of Eqs.(C-18) and (C-22). 

Let us exploit in the second step the evolution equations for mode functions to generate 
solutions with a finite value of the three- momentum q in addition to uo. This corresponds 
to using the mode function w^ k (x, x ) = e* q ' XT " ! '-u -k( x , #o)- We obtain the following general 
expression : 



S(p) = (2n( Po ) + 1) 



(po - UT y ) 2 - (p - qr y ) 2 - m? - ^r z + ie 



(D-14) 



We can get .M? 2 (x, x) or A^^(k), which are introduced by .Mi 2 (x, x ) = Sk-Mi 2 (k), per- 
forming the inverse Fourier transformation from S(p) by the use of the relation of Eq.(D-6) 



1 



(D-15) 



(k - UTy) 2 - (k - qr y ) 2 - m 2 - \j^t z + ie ' 

Let us consider the propagator in the subspace of the first two isospins. First, we con- 
centrate on the time-like q 2 as q 2 = u 2 with q = 0. After integration with respect to ho, we 
have Mf 2 (k) for a,b = 1,2 : 



(£i(x)0i(x)) = .M^(x,x) = 5 n (x,x) 

=/^W )+ i) 



1 + 2^ 2 + /i2 



+ (2n(£_) + 1) 



4E+ 2E + (E\ — E 2 ) 
1 2u 2 + /ig 



(<£ 2 (x)<£ 2 (x)) = .Mi 2 (x,x) = <5 22 (a;,a;) 
d 3 k 



/ 



(27T) 



(^ a (x)^ ft (x)) = 



(2n(£+) + 1) 

+ (2n(£_) + 1) 
for a 7^ b , 



2E_{E\-E 2 _) 

2u 2 - //§ 
2E + (E 2 + -E 2 _ 

2u 2 - $ 



AE 4 



+ 



4£_ 2E_{E\-E 2 l) 



where we denote (x, x^S\y, y ) as S(x,y). Here, E ± are defined in (C-7) and 

1 



n(E) 



(D-16) 



(D-17) 



e^ - 1 ' 

Also, (y3 c (x)y3 c (x)) for c > 3 are given in Eq.(D-2). Thus, our assumption that M 2 b has a 
diagonal form is valid because (</? a (x)</? 6 (x)) = for a 7^ 5. 
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Secondly, let us consider the space-like q 2 as q 2 = — q 2 with u = 0. In this case, we can 
get A^i2( x ; x ) m the same way as the time-like case : 



(£i(x)0i(x)) = Mn(x,x) = S n (x,x) 



+ 



(4 



+ (2n(E_) + 1) 



4E + 2E + (E 2 + -E 2 _) 
1 



Hi 



4E_ 2E_{E\ — E 2 ) 



(^ 2 (x)^ 2 (x)) = .M 22 (x,x) = 522(0:, X) 

<J ''' k '(2„(A+)4 1) 



(27T) 



(^o(x)^(x)) = 



+ (2n(£_) + 1) 
for a 7^ 6 . 



4£+ 2E + (E\-E 2 _\ 
1 



+ 



4£_ 2E_{E\-E 2 _) 



Here, (> 0) are obtained by the poles of the integrand in (D-15) with uj = : 



E 2 ± (k) = k 2 + q 2 + M 2 ± J ^ + 4(k • q 



(D-18) 



(D-19) 



Finally, the gap equation in the isospin-rotating frame at finite temperature is derived 
from the equation of motion for (p with (p a = (poS a i. We obtain the gap equation for y? in 
the symmetry broken phase : 



A A N A 

-ifl = uj 2 - q 2 - m 2 - - Y^{0 c (x)0 c (x)) - -(^i(x)^i(x)) 



(D-20) 



due to (</? a (x)</? 6 (x)) = for a 7^ 6. By using the mass M 2 , the above gap equation is recast 
into simpler form as 

(D-21) 



-ifl = -uj 2 + q 2 + Ml 
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